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, Abstract. It has been shown by H. Bethe more than 70 years ago that the number of excited states of a Fermi gas grows, 

O |. at high excitation energies Q, like the exponential of the square root of Q. This result takes into account only the average 

density of single particle (SP) levels near the Fermi energy. It ignores two important effects, namely the discreteness of the 
SP spectrum, and its fluctuations. We show that the discreteness of the SP spectrum gives rise to smooth finite-<2 corrections. 
Mathematically, these corrections are associated to the problem of partitions of an integer. On top of the smooth growth of 
the many-body density of states there are, generically, oscillations. An explicit expression of these oscillations is given. Their 
properties strongly depend on the regular or chaotic nature of the SP motion. In particular, we analyze their typical size, 
temperature dependence and probability distribution, with emphasis on their universal aspects. 
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The main theoretical framework to understand the behavior of the density of states of a Fermionic system has been, 
and still is, the independent particle model. In the degenerate gas approximation, when the excitation energy of the gas 
is much smaller than the Fermi energy £ F , the excitation spectrum relies on the properties of the single-particle (SP) 
spectrum near s F . In most theoretical calculations only the average SP density of states, p, is taken into account. For a 
system of A noninteracting fermions moving in a mean-field potential, the number of excited states of the many-body 
(MB) system contained in a small energy window dQ at energy Q, PMB{A,Q)dQ, is 1 1 j 
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INTRODUCTION 



pM S (A,e) = ^=L-exp(2 v / fl e) , a = 7t z p/6. (1) 
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Here Q is measured with respect to the ground state energy of the gas. For a given potential, p is in general a function 
of A. Besides the condition Q <C E P , this expression assumes that the excitation energy is large compared to the SP 
mean level spacing 8 = p ~ 1 . For simplicity, and because our aim here is not to compare with experimental data, we 
ignore angular momentum conservation, isospin, etc. 

After Bethe's work, more accurate calculations of p^g were made 0]. Schematic shell corrections related to a 
periodic fluctuation of the SP density were computed in 1 3] (see also |4j), introducing the so-called back-shifted Bethe 
formula. Several phenomenological modifications of Eq.Q have been proposed to match the experimental results. 
These models take into account, for example, shell effects, pairing corrections and residual interactions 01 0,0,01, 
introducing a multitude of coexisting phenomenological parameterizations. The present status of the understanding do 
not allow to draw a clear theoretical picture of the functional dependence of p MB with A and Q. 

Our purpose is, within a SP picture, to further develop the theoretical analysis to include several important effects 
that are missing in Eq.([0. This presentation is based on the work reported in |9[. 

There are different ways in which Eq.Q can be improved. As mentioned before, it is the leading term of an 
expansion valid for a large number of particles A and for high energies Q (compared to the SP mean level spacing 
8). For instance, assuming the gas is confined by a Woods-Saxon like mean field potential, then 

3 A 7t 2 A , 

p = --^ fl = -A,-MeV->, (2) 
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FIGURE 1. The parameter a appearing in the Fermi gas level density formula as a function of the mass number A (taken from 
Ref.H). 

where we used the approximation e F « 37 MeV. This expression of p is the leading order term of an expansion in 
decreasing powers of A. Corrections to it can be incorporated by considering lower order terms in the Weyl series, 
that includes surface corrections, curvature corrections, etc 1 10]. A similar expansion in decreasing powers of A holds 
for the energy or mass of the nucleus (the liquid drop formula). As is well known, the coefficient of each term in the 
expansion of the energy obtained from a Fermi gas is not correct, since the renormalization due to interactions is out 
of the range of the model. Similarly, the coefficient 1/15 in Eq.(0 does not give a good description of the average 
trend of the experimental data obtained from neutron resonances, which is closer to 1 /8 (cf Fig. 1 ). Although of great 
interest, we will not discuss these A corrections and their renormalization due to interactions, but concentrate on two 
other corrections that are within the scope of a noninteracting model. The first one (next section) leads to smooth 
lower order terms in the excitation energy Q. The second type of corrections incorporates oscillatory terms in A and 
Q. These oscillations, superimposed to the smooth A/8 trend mentioned above, are clearly visible in the experimental 
data shown in Fig.l. 



DISCRETENESS OF THE SPECTRUM: AVERAGE BEHAVIOR 

The only property of the SP spectrum in Eq.Q is the average density of levels p. This approximation treats the SP 
spectrum as a continuum, and ignores the influence in the many-body density of states of the discreteness of the 
spectrum and of the exact position of the SP levels. One step beyond, that incorporates the discreteness but ignores the 
fluctuations, is to consider a locally perfectly regular SP spectrum (ID harmonic oscillator), given by the SP energies 
£„ = n8, where n is an arbitrary integer, and 8 is the distance between neighboring levels. The Fermi gas model consists 
of A noninteracting fermions that occupy the equidistant levels, with occupation number or 1 . In the ground state 
the particles fill the lowest A SP states. We restrict the analysis to the approximation Q <C e F , in which the number of 
excited particles is small compared to the total number A, and effects due to a finite number of particles can be ignored. 
In the excited states particles occupy SP levels above e F , thus creating holes in the Fermi sea. The possible excited 
energies, measured with respect to the ground state energy, are Q = m8, where m is an arbitrary positive integer. The 
excitation energies are thus trivial. The nontrivial information comes from the degeneracy of each of theses energies, 
since there are many different many-body configurations with the same excitation energy. The problem then is how to 
compute the degeneracy of each excited state of energy Q — m.8. 

The solution is obtained by realizing that the problem we are facing is exactly equivalent to a well known problem 
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FIGURE 2. Ground state (G.S., left column) and excited states of energy Q/& =4 (next five columns) of a gas of noninteracting 
fermions in a perfectly regular SP spectrum. Each excited state is identify with a decomposition of the excitation energy 4 as a sum 
of integers (bottom of each column). 



in mathematics, i.e. in how many different ways an integer can be decomposed as a sum of integers (partitions of an 
integer, p(m)). We will simply provide here a graphical illustration of this equivalence (Fig. 2). Consider for example 
the case m = 4, the argument easily generalizes to any m. m = 4 can be decomposed as a sum of integers in 5 different 
ways (and thus p(4) =5): 1 + 1 + 1 + 1,2 + 1 + 1,2 + 2, 3 + 1 and 4. As Fig. 2 shows, each of these decompositions 
is in a one-to-one correspondence with one of the possible excited states of the Fermi gas of total energy Q/8 =4. 
The ground state is shown on the left column. In the second column from the left four particles have been pushed one 
level up, hence the equivalence with 1 + 1 + 1 + 1 (other mappings may be used as well). In the third column the upper 
particle was pushed up by two, the next two by one. Etc. As has been emphasized many times, there is also a (trivial) 
one-to-one correspondence between the partitions of an integer and the excited states of a gas of noninteracting bosons 
in a one-dimensional harmonic oscillator potential. In this case each integer in the decomposition of the excitation 
energy indicates the level occupied by one of the excited particles. 

The partition of an integer grows very fast with m: it is 5 for m — 4, whereas p(2QQ) = 397999029388 (this number 
was computed by hand by P. MacMahon in 1918). The generating function of p(m) was obtained in 1753 by L. Euler. 
An explicit expression for p(m) came much later, under the form of an asymptotic (exact) formula (initially obtained 
by Hardy and Ramanujan 1 1 1 ], improved and made rigorous later on by Rademacher 1 12]). To facilitate the comparison 
with more general approaches, we express the result in terms of ~pQ, where here p = 5 _1 is the inverse of the exact 
spacing between neighboring levels. The first terms of the expansion are given by, 

log {pSS/p) =^\% 2 PQ~ log ( V48 pQ) - ^^(p fi )-i_^_^(p Q )-i + ^((p Q )-!). (3) 

We emphasize with the notation p^g that this is the many-body density of states for a Fermi gas in a ID harmonic 
oscillator spectrum. As expected, the first two terms of this equation reproduce Bethe's formula. However, Eq.(|3} goes 
much further: it is a truncation of the exact asymptotic expansion of the density of states, whose precision can be 
improved by adding further terms. Notice that all the terms in Eq.Q are smooth functions of the excitation energy Q. 

It is natural to ask about the relevance of Eq.0 in realistic systems, where the true SP spectrum consists of discrete 
energy levels arranged, in most cases, with no particular order. Equation Q clearly goes beyond Bethe's result, since 
it describes the exact behavior of the many-body density of states of a gas of particles that occupy a discrete, perfectly 
regular arrangement of SP levels. It is therefore not unreasonable to expect that the result is valid for an arbitrary 
spectrum having the same average density p, and that the effects of the fluctuations of the SP energy levels with 
respect to a perfectly regular arrangement come on top of it. Although we do not have for the moment an explicit proof 
of this statement, numerical simulations to be presented below seem to confirm this hypothesis. 



FLUCTUATIONS 



The exact dependence of p MB on A and Q is sensitive to the detailed arrangement of the SP energy levels around the 
Fermi energy. Strong deviations with respect to a regular spacing, with possible degeneracies, may be produced by the 
presence of symmetries of the confining self-consistent potential. These deviations induce, in turn, oscillations in the 
thermodynamic functions of the gas. Shell effects are therefore due to deviations (or bunching) of the SP levels with 
respect to a perfectly regular spectrum. The degeneracies of the electronic levels of an atom produced by the rotational 
symmetry are an extreme manifestation of level bunching. In general, in systems that have other symmetries, or have 
no symmetries at all, there will still be level bunching, but its importance will typically be minor. Therefore, depending 
on the presence or absence of symmetries, the shell effects may be more or less important. The level bunching, and 
more generally the fluctuations of the SP energy levels, are thus a very general phenomenon. The theories that describe 
those fluctuations make a neat distinction between systems with different underlying classical dynamics (e.g., regular 
or chaotic). Our understanding of these fluctuations and of their connections with the regular or chaotic nature of the 
SP motion has greatly increased in the last decades. For a recent review see 0. 

Oscillations in the many-body density of states of nuclei are clearly visible in the experimental data available from 
neutron resonances as a function of the mass number A. The parameter a extracted from the data, plotted in Fig.l, 
shows oscillations around an average growth (A/8). The effective description of the oscillations through the parameter 
a is an artifact of the analysis, since the Fermi gas model relates a to p , a smooth function. As we will see below, the 
theoretical analysis leads to oscillatory corrections that enter as an additional term in the exponent of p MB , and thus 
cannot in general be interpreted as an effective a. 

There are two distinct and important scales to describe the SP fluctuations. The first one, the smallest one, is the 
mean spacing between SP energy levels 8. The fluctuations of the SP energy levels on that scale have been shown 
to be, for a large class of systems, universal. Their statistical properties are described by uncorrelated sequences for 
integrable systems, and by random matrix theory in chaotic ones. The second relevant energy scale is E c — /i/T m i n , 
where is of the order of the time of flight across the system (a precise definition will be given below). For A large, 
the ratio g = E c / 8 is much larger than one (for instance, g ~ A 2 / 3 in a three dimensional cavity). In contrast to the 
fluctuations on scales 8, the fluctuations on scales E c are long range modulations of the SP spectrum whose structure 
and amplitude are system specific. There is therefore no universality on this scale. 

For an arbitrary SP spectrum the computation of the density of excited states of a fermionic system is a difficult 
combinatorial problem for which no exact solution is known. What we are seeking here are the variations of the MB 
density of states due to fluctuations of the SP spectrum with respect to a perfectly ordered spectrum, described in the 
previous section. An explicit formula that includes these effects has been obtained recently |9]. The result, obtained 
from a saddle point approximation of an inverse Laplace transform of the MB density, takes the form |9] 



Pmb(A,Q) = 



1 
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exp 




(4) 



The function §(\i , T) is the fluctuating part of the energy of the gas at chemical potential ji and temperature T, defined 
as follows. The energy is, as usual, given by 



&(H,T) = dee p(e) 



1 



(5) 



where the function p(e) is the SP density of states defined by the SP energies £; 



The SP density of states is usually decomposed into a part p that has a smooth dependence on the energy e plus another 
contribution p that describes the deviations with respect to the smooth behavior. The fluctuating part of the energy is 
defined as 



£{li,T) = g(li,T)- / deep{e 



l +e (.e-n)/T 



(6) 



In Eq. (0} the arguments of § are functions of A and Q, ji = ji(A,Q) and T = T(A,Q). These two functions are 
determined by inversion of the equations 



A = JT{H,T), 

Q = ^(M,7-)-^(M,0). 



(7) 



The function JY{n,T) = Jds p(e) 



i-i 

1 +e( £ -V)/ T i s the particle number. We are ignoring here (small) variations 

of the chemical potential with temperature. To lowest order, in which only the average behavior of the energy and of 
the particle number is taken into account, the second equation in Q leads to the usual relation between temperature 
and excitation energy 

' Q (8) 

and pi is simply a function of A, jx = ju(A). For instance, for a Fermi gas in a 3D cavity of volume V and particle mass 
m, 

2nfi 2 /3V7rA\ 2/3 



^^r\r^r) ■ (9) 

Equation @ shows that, superimposed to the smooth growth of the density described by Eq.Q, there are oscillatory 
corrections in the exponent of the density of states. These corrections are directly related to energy fluctuations of the 
system. 

The oscillatory nature of the new corrections is clearly displayed through a semiclassical theory, that expresses the 
quantum properties of the fermion gas in terms of classical solutions of the SP equations of motion. In this approach 
S"(jj.,T) is written as a sum over the classical periodic orbits of the mean-field potential 013 El 



r) = in 2 V V Ap > r k }:* p/Tt) cos (r s p /h + v p , r ) . 



P r=l r l P 



(10) 



Each orbit is characterized by its action S„, stability amplitude A p/ , and Maslov index v p . r . t p is the period of the 
periodic orbit, and k{x) —xj sinh(x) is a temperature factor that introduces the time scale ty = h/(nT) conjugate to the 
temperature. This expression describes the departures of S with respect to its mean behavior due to the fluctuations of 
the SP spectrum. The behavior of § strongly depends on whether fiorT (or A or Q) are varied. A temperature variation 
modifies the prefactors of the summands in S (through the function k), and therefore produces gentle variations of the 
fluctuating part in the exponent of the density of states 



s(a,q) = <'/(,u.os-/;(,u.r: 



(ii) 



In contrast, A p ., -, x p and S p depend on [i (and therefore on A). For large values of fi, S p ^ Ti and the dominant variation 
with the particle number (or any other parameter that modifies the actions) comes from the argument of the cosine 
function in £ . Rapid oscillations of S are therefore generically expected when the number of particles is varied. 

For a given mean field potential there exists an infinite number of periodic orbits p. The spectrum of periods x p has 
no upper bound. It has, however, a lower bound, given by the period T m i n of the shortest periodic orbit. Since T m i n is 
usually the smallest characteristic time scale in the system, it determines the largest energy scale E c — hjXmm in which 
modulations (bunching) of the SP energy levels occur 

If the classical periodic orbits are known in a particular system, the corrections can be computed explicitly. The 
periodic-orbit sum is dominated by the short orbits, and rapidly gives a good approximation of the result. To illustrate 
our findings we have tested some of the predictions by a direct numerical counting of the MB density of states in a 
particular system. Figure 3 shows the results obtained for a gas of about 2000 fermions contained in a two-dimensional 
rectangular cavity, an integrable system. For each number of particles we compute the MB density of states at three 
different temperatures, measured in units of 

T s = 8/2n 2 and T c =E c /2% 2 . 

The theoretical curve is computed according to the expression 

log(p MB /p) = log(p™/p) + S(A,Q) , (12) 



where S is calculated from the periodic orbits of the rectangle. 

If the smooth corrections to the density of Eq.Q are not included, and only the result of Bethe is used, a systematic 
deviation is observed in the average behavior of p MB between theory and numerics. Although the theoretical analysis 
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FIGURE 3. The logarithm of the many-body density of states of A noninteracting fermions in a rectangular billiard of sides 



a = y (1 + v / 5)/2 and = a 1 . The value of g is 100. Dots: numerical computation at three different excitation energies; solid 
curves: theoretical prediction 1 1 21 : dashed curves: smooth part Eq.J5J. 



of Ref. 1 9] does not include those corrections, the numerical results seem to indicate that their validity goes beyond the 
simple ID harmonic oscillator spectrum. Notice the extremely good accuracy of Ea.l fT2l . either for the average value 
of the density as well as for the fluctuations. Notice also that the typical size of the fluctuations can be quite important 
(it is of the order of the average for the lowest excitation energy computed). 



STATISTICAL ANALYSIS OF THE FLUCTUATIONS: UNIVERSALITY 

Equation H2i decomposes the exponent of the density of states into a smooth part (given by the Hardy-Ramanujan 
result) plus oscillatory contribution, S(A, Q), associated in a semiclassical theory to a sum over periodic orbits. S(A,Q) 
relates the shell corrections in the MB density of states to the fluctuations of the energy of the gas. § depends on 
temperature only through the function k (cf Eq.fTol'). The effect of this function is to exponentially suppress the 
contribution of periodic orbits whose period z p T T 1161 1 1 511 - Since there is no suppression at T = (because 
x T — > °o), only orbits whose period T p ~ T/ r contribute to the difference ^(/i,0) — ${ll,T). At temperatures such that 

%j <C T m i n , the term £hi,T) becomes exponentially small, and only <f(/i,0) remains. The shell correction S(A,Q) 
therefore decays as T~ l atT ^> T r . This decay contrasts with the more common exponential damping observed in 
other thermodynamic quantities Hal . 

A clearer picture of how the fluctuations behave may be obtained through a statistical analysis. The most simple 
property of the fluctuations is (5) = 0, where the brackets denote an average over a suitable chemical potential (or 
particle number) window. This result is valid only to first order in the expansion leading to Eq.@ |9]; it can be shown 
that higher order terms contribute to a non-zero average. The next non-trivial statistical property is the variance (S 2 ), 
that may be computed using Ea. (llO> . The result is (S 2 ) = (1 /2) JqcIx K(x,xn) [1 — fc(x)] 2 /x 4 , where K(x,xn) is the 
rescaled form factor of the SP spectrum (cf Eq.(36) in Ref.| 1511 ). The latter function depends on therescaled Heisenberg 
time Xft = hp '/z T . It describes system-dependent features for x of the order of x m ; n = T T , while it is believed to 
be universal forx ^> x m ; n . The universality class depends on the regular or chaotic nature of the dynamics, and on its 
symmetry properties. 




Taking into account the basic properties of K(x,Xn), in chaotic systems three different regimes for (S 2 ) as a function 
of temperature are found |9] (remember the relation © between T and Q): 

• Low temperatures T <C Tg, In this regime 

(S 2 )=c 4 T/(2T s ) , (13) 
where C4 = 0.0609... is the value at k = 4 of the constant 



Ck = 



x 



d x 



1- 



sinhx 



(14) 



• Intermediate temperatures Tg <c T <C T t . In this regime the size of the fluctuations saturates at a universal constant 

(S z )=c 3 /p, (15) 

where C3 = 0. 1023 . . . and j8 = 1 (2) for systems with (without) time-reversal invariance. 

• High temperatures T ^>T C . The size of the fluctuations decreases with excitation energy, 

(5 2 ) = <J%,0))[l-8e- r ^]/r 2 . (16) 

After an exponential transient, a power-law decay (S 2 ) 1 ' 2 « T~ l is thus obtained. 

The situation is different in integrable systems, where only two regimes are found. At low temperatures the result 
is identical to that of chaotic systems. The linear growth at low temperatures is thus totally universal and independent 
of the system. The difference is that in integrable systems the growth extends up to much higher temperatures, 
T « T c , without saturation. At that temperature the variance of the fluctuations is of order g. In integrable systems, 
the maximum amplitude of the fluctuations is therefore reached at T « T c , and its typical size is much larger than in 
chaotic systems. At high temperatures T ^> T c the decay is almost identical to that of chaotic systems (the coefficient 
8 is replaced by a 12). A schematic representation of the temperature dependence of the variance of the fluctuations 
for chaotic and integrable systems is given in Fig. 4. 

It has been shown 11511 that «?(ju,r) is dominated, at any T, by the shortest classical periodic orbits. In contrast, 
the difference <?(ju,0) — <?(ju, T) depends on orbits whose period % p ~ % T . For temperatures T -C T c the statistical 
properties of these orbits are universal, and correspondingly the probability distribution function of S is expected to 
be universal, in the sense that at a given temperature it should only depend on the nature of the underlying classical 
dynamics (regular or chaotic), and the symmetries of the system. This statement is supported by the fact that in the 
limit T — > 0, S w —d^dx^/dT. The probability distribution of the latter quantity was studied in Ref. lfl5Tl : it was 
shown that it coincides at low temperatures with that obtained from a Poisson spectrum for integrable systems and from 
a random matrix spectrum for chaotic ones. This confirms the universality of the distribution at low temperatures. As 



the temperature is raised, the universality of the probability distribution of S will be lost for temperatures of the order 
or greater than E c , where system specific features (i.e., short periodic orbits) are revealed. In this respect, the decay 
(I16> is only indicative. Its exact form depends on the details of the short periodic orbits spectrum, treated here only 
through a rough approximation. 



CONCLUDING REMARKS 

Two improvements with respect to Bethe's formula of p MB have been discussed. The first one incorporates the 
discreteness of the SP spectrum by considering a set of equidistant SP states (ID harmonic oscillator). The exact 
many-body density of states is obtained by mapping the problem to the computation of the number of decompositions 
of an integer as a sum of integers. The exact formula adds finite-2 corrections to the asymptotic result. Although a 
proof is missing, numerical simulations suggest that the range of validity of this result extends beyond the ID harmonic 
oscillator spectrum to any spectrum, in the sense that it describes with good accuracy the smooth part of the growth of 
Pmb in systems with an arbitrary SP spectrum (cf Fig. 3). 

Beyond the smooth behavior, the second improvement concerns the fluctuations or shell effects in p MB . In a 
semiclassical theory, each periodic orbit has been shown to contribute to log (Pms/p) with a fluctuating term as a 
function of the number of particles A, of wavelength 

AA = ^. 
8 

At low excitation energies long orbits contribute to the fluctuations while short ones are exponentially suppressed, 
leading to wild oscillations and universality of their statistical properties. As the temperature increases, the situation 
is reversed, long orbits are exponentially suppressed while short orbits come into play. At T = T c , oscillations of 
wavelength of order AA = E c /8 — g are predicted. These oscillations are clearly visible, with the correct wavelength, 
in Fig. 3. With the raise of the short periodic orbits as T increases, the universality of the statistical properties of the 
fluctuations disappears for T of the order or higher than T c . 

Concerning the typical size of the oscillations, a nontrivial dependence as a function of excitation energy (or 
temperature) was found. The results are summarized in Fig. 4. The behavior is quite different in chaotic and integrable 
systems. At low temperatures the variance grows linearly with T in all cases. However, a plateau is rapidly reached at 
T = Tg in chaotic systems, whereas the amplitude of the fluctuations continue to grow in regular ones. At T = T c the 
differences reaches its maximum amplitude, where the variance in regular systems is g times larger than in chaotic ones 
(remember, moreover, that g ~ A 2 / 3 in a three dimensional cavity, so that the difference increases with an increasing 
number of particles). At this point the situation is similar to what was found for the ratio of the variance of the 
fluctuations of the nuclear mass due to regular and chaotic motion ifTvtl . At T ^ T c the variance decreases as T~ 2 for 
both types of dynamics. The validity of the results discussed here in the analysis and interpretation of the experimental 
data on the nuclear level density will be presented elsewhere. 
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